Phase transitions of quasistationary states in the Hamiltonian Mean Field model 
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The out-of equilibrium dynamics of the Hamiltonian Mean Field (HMF) model is studied in 
presence of an externally imposed magnetic field h. Lynden-Bell's theory of violent relaxation is 
revisited and shown to adequately capture the system dynamics, as revealed by direct Vlasov based 
numerical simulations in the limit of vanishing field. This includes the existence of an out-of- 
equilibrium phase transition separating magnetized and non magnetized phases. We also monitor 
the fluctuations in time of the magnetization, which allows us to elaborate on the choice of the 
correct order parameter when challenging the performance of Lynden-Bell's theory. The presence 
of the field h removes the phase transition, as it happens at equilibrium. Moreover, regions with 
negative susceptibility are numerically found to occur, in agreement with the predictions of the 
theory. 



PACS numbers: 

Long-range interacting systems are characterized by a 
slowly decaying interparticle potential, which in fact re- 
sults in a substantial degree of coupling among far away 
components. In these systems, energy is consequently 
non-additive and this fact yields a large gallery of pe- 
culiar, apparently unintuitive, phenomena: the specific 
heat can be negative in the microcanonical ensemble, 
and temperature jumps may appear at microcanonical 
first-order phase transitions [ll, Canonical and mi- 
crocanonical statistical ensembles can therefore be non- 
equivalent in presence of long-range interactions, an in- 
triguing possibility which has been thoroughly discussed 
working within simplified toy models. 

Systems subject to long range couplings also dis- 
play unexpected dynamical features. Starting from out- 
of-equilibrium initial conditions they are occasionally 
trapped in long lasting regimes, termed Quasi Stationary 
States (QSS), whose lifetime diverges with the number of 
elements, N, belonging to the system under scrutiny Q. 

The QSSs have been shown to relate to the stable 
steady states of the Vlasov equation, which governs the 
dynamical evolution of the single particle distribution 
function in the continuum limit — > cxo [J [^-[^ . Work- 
ing within this setting, one can implement an analytical 
procedure, fully justified from first principles, to clarify 
some aspects of QSS emergence. The idea, inspired to the 
seminal work of Lynden-Bell is based on the defini- 
tion of a locally-averaged ( "coarse-grained" ) distribution 
function, which translates into an entropy functional, as 
follows from standard statistical mechanics prescriptions. 
By maximizing such an entropy, while imposing the con- 
straints of the dynamics, returns a closed analytical ex- 
pression for the single particle distribution function of 
the system in its QSS regime. The predictive adequacy 
of the technique was tested versus numerical simulations 



for specific applications relevant in e.g. plasma physics, 
and for the Hamiltonian Mean Field (HMF) model Q , to 
which we will make extensive reference in the following. 
Furthermore, the Lynden-Bell approach allows one to 
successfully identify out-of-equilibrium phase transitions 
separating homogeneous and non homogeneous steady 
states [1, @- More recently, the Lynden-Bell procedure 
was applied to an open HMF system, modified by the 
inclusion of an externally imposed field, prognosticating 
the existence of regions of negative susceptibility which 
were then observed in direct simulations of the discrete 
Hamiltonian [Toj . 

In this paper we revisit the Lynden-Bell analysis for 
the HMF model. The theoretical scenario is tested ver- 
sus Vlasov based simulations, returning an overall good 
agreement. We discuss also the impact of the choice of a 
monitored quantity on the characterization of the order 
of the phase transition in absence of the external field. 
The role of the externally imposed field is assessed, with 
emphasis on the modification of the Lynden-Bell transi- 
tion. The response of the system to the external forcing 
results in a smoothing of the transition that separates 
homogeneous and non homogeneous regimes, an obser- 
vation which a posteriori supports the identification of 
such phenomenon with a genuine phase transition. We 
here anticipate that regions with negative susceptibility 
will be also identified in agreement with the Lynden-Bell 
scenario depicted in [loj . 

The paper is organized as follows. The next section is 
devoted to introducing the HMF model and to discussing 
its continuous analogue. In Section |ll] we present the 
Lynden-Bell calculation, revisiting the results with refer- 
ence to the unforced system. In Section HIIl we present the 
results of the numerical simulations, based on a Vlasov 
code, aimed at verifying Lynden-Bell's prediction of the 
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presence of an out-of-cquilibrium phase transition in the 
HMF modeL The effect of applying an external magnetic 
field h is discussed in Section |lVl Finally, in Section Ivl 
we sum up and conclude. 



I. THE HAMILTONIAN MEAN FIELD MODEL 

The Hamihonian Mean-Field (HMF) model de- 
scribes the motion of N classical rotors coupled through a 
mean-field interaction. The system, in its standard for- 
mulation, can be straightforwardly modified to include 
an external perturbation that acts on the particles as 
a magnetic- like field [ll|. The model is mathematically 
defined by the following Hamiltonian: 

N N N 

^=2^^^'+ 27V ^ [l-cos{9,~e,)]^hY,co8{e,), 
" i=i ij=i j=i 

(1) 

where 6j represents the orientation of the j-th rotor and 
Pj is its angular momentum. The scalar parameter h 
measures the strength of the magnetic field. Hamiltonian 
dl]) with h = has been widely studied in the past as a 
prototype model of long-range interacting systems. To 
monitor the dynamics of the systems, one often refers to 
the magnetization, a collective variable defined as 

1 ^ 

m= —'^{cos 9i, sin 9i) ^ {rrij;, my) . (2) 
1=1 



The modulus of m, m ~ + m^, measures the de- 

gree of bunching of the rotors along a given direction. 
The model can be also interpreted as describing iV parti- 
cles moving on a circle. Within this interpretation, mag- 
netized regimes signal the presence of a localized agglom- 
eration of particles on the circle. 

As previously reported in the literature Q, starting 
from an out-of-equilibrium initial condition, the system 
gets trapped in long lasting QSSs, whose macroscopic 
characteristics differ significantly from those associated 
to the corresponding equilibrium configurations. QSSs 
develop for both /i = and h settings. In the N ^ oo 
limit, the system is indefinitely stuck in the QSS phase. 

On the other hand, when performing the limit for infi- 
nite system size, the discrete model Hamiltonian ([T|) ad- 
mits a rigorous continuous analogue. This is the Vlasov 
equation which governs the evolution of the single parti- 
cle distribution function f{9,p,t): 

dl df dv[f] df ^ 

dt ^86 dO dp 
V[f]{6) = I ~{m,[f] + h) cose - my\f]sine , 
[/] = J dOdpfcosO , 

rriyif] = [ de dp f sine, (3) 



where V is the interaction potential that depends self- 
consistently on f{e,p,t). According to this kinetic pic- 
ture the free streaming of the particles is opposed by 
a potential term V[f], reminiscent of the discrete for- 
mulation, expressed clS &- self-consistent function of the 
dynamically varying distribution f{e,p,t). 

In light of the above, the QSSs have been interpreted 
as stable steady states of the underlying Vlasov equation. 
Working within this setting, and invoking the aforemen- 
tioned Lynden-Bell violent relaxation theory @ , one can 
progress analytically at least for a simplified choice of the 
initial condition. A short account of the technicalities is 
provided in the following section. 



II. THE MAXIMUM ENTROPY SOLUTION 

Assume the particles to be confined within a bounded 
domain of phase space, therein displaying a uniform 
probability distribution. Label /o the constant value 
of f(e,p,t) within the selected domain, as imposed by 
the normalization condition. This working ansatz corre- 
sponds to dealing with the "waterbag" distribution: 

f /o = iski if - Ap < p < Ap 

f{e,p, 0) = < &nd -Ae <e < Ae (4) 

[ otherwise, 

that is even in both e and p: f{—e,—p,0) = f{e,p,0). 
For distributions endowed with this symmetry, it can 
be straightforwardly proven that, being m,y = ini- 
tially, its value remains zero during time evolution. This 
in turn implies that also the total momentum P = 
J pf{e,p, t)dedp, which is zero initially, remains zero dur- 
ing the whole time evolution, i.e. there is no global ro- 
tation of the particles on the circle. With this choice, 
one parametrizes the initial condition in terms of the 
energy density u = H/N and the initial magnetization 
m = (77io,0). Momentum P cannot be considered as a 
global invariant, because the presence of an external mag- 
netic field breaks the translation symmetry e ^ e + a. 
However, for the initial distributions Q, momentum is 
fixed to zero. 

Under the Vlasov evolution, the waterbag gets dis- 
torted and filamented at smaller scales, while preserv- 
ing its surface in phase space. The distribution stays 
two-level (0, /o) as time progresses. By performing a 
local average of / inside a given mesoscopic box, one 
gets a coarse-grained profile which is hereafter labelled 
/. As opposed to /, the locally averaged function / con- 
verges to an asymptotic equilibrium solution which can 
be explicitly evaluated via a rigorous statistical mechan- 
ics procedure, adapted from the pioneering analysis of 
Lyndcn-BcU. An entropy functional s(/), can be in fact 
associated to /, through a direct combinatorial calcula- 
tion 0. In the two-level waterbag scenario, the mixing 



entropy density s takes the form: 
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The energy density 
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is conserved. In addition, the normahzation of the distri- 
bution / has to be imposed, which physically corresponds 
to assume constant mass. Requiring the entropy to be 
stationary, while imposing the conservation of energy and 
mass, result in a variational problem that admits the fol- 
lowing solution: 



fo 



I -j- g/3./o(p"/2 — rrij. cos 9 — iriy svnO — h cos 9) + a 



(7) 



where a and /? are Lagrange multipliers associated, re- 
spectively, to mass and energy conservation and mx and 
niy depend on /qss- The self-consistent nature of Eq. ([7]) 
is evident: mx[f qss] and mj^[/Qss] are functional of /qss 
and both enter in the determination of /qss itself. We 
also emphasize that the label QSS is introduced to re- 
call that the stationary solution of the Vlasov equation 
are indeed associated to QSSs of the discrete A^-body 
dynamics. 

Back to solution ([7]), one can determine the predicted 
values of mx, my, a and (3 once the energy e, the field 
h and the waterbag height fo are being assigned. This 
step is performed numerically, at sought accuracy, via a 
Newton-Raphson method. 

Consider first the limiting case h = 0. Depending on 
the value of the predicted magnetization m, one can ide- 
ally identify two different regimes: the homogeneous case 
corresponds to to = (non- magnetized), while the non- 
homogeneous setting is found for (magnetized) m ^ 
solutions. A phase transition 0, [H materializes in the 
parameters plane (mo,it) and the resulting scenario is 
depicted in Fig. [1] When fixing the initial magnetization 
and decreasing the energy density, the system passes from 
homogeneous to non-homogeneous QSS. The parameters 
plane can be then formally partitioned into two zones 
respectively associated to an ordered non-homogeneous 
phase, m ^ 0, (lower part of Fig. [J), and a disordered 
homogeneous state, m = (upper part). These regions 
are delimited by a transition line, collection of all the 
critical points (to^,?!*^), which can be in turn segmented 
into two distinct parts. 

The full line stands for a first order phase transition: 
the magnetization experiences a finite jump when cross- 
ing the critical value (mg,?/^). Conversely, the dashed 
line refers to a second order phase transition: the magne- 
tization is continuously modulated, from zero to positive 
values, when passing the curve from top to bottom. First 
and second lines merge in a tricritical point. 



iL^ L 



0.04 0.06 0.08 1 0.12 014 0.16 




1 order 

2 order 

% tricritical point 



FIG. 1: Phase diagram in the {mo,u) plane at /i = 0. The full 
line refers to the the first order transition, while the dashed 
line stands for the second order one. The symbol traces the 
exact location of the tricritical point. Inset: lateral edges of 
the coexistence regions in the first order region, as predicted 
by the theory. 



The case with h ^ has been recently addressed in 
[loj for what concerns the Lynden-Bell theory and work- 
ing at constant fo, while th e eq uilibrium properties have 
been thoroughly studied in [11|. Again, the Lynden Bell 
approach proves accurate in predicting the macroscopic 
behavior as seen in the iV-body simulations. Interest- 
ingly, below a threshold in energy the system shows 
negative susceptibility x = dM/dh, the magnetization 
decreasing when the strength of h is enhanced. Con- 
versely, above the critical energy value, the magnetiza- 
tion amount grows with h, which corresponds to deal- 
ing with positive susceptibility. Besides providing an a 
posteriori evidence on the adequacy of the Lynden-Bell 
technique, the presence of a region with negative sus- 
ceptibility, was interpreted in (lO| as the signature of 
an out-of-equilibrium ensemble inequivalence. Further- 
more, the presence of the field h removes the phase tran- 
sition and the magnetization continuously decreases from 
unity, at zero temperature, to zero, at infinite tempera- 
ture. Therefore, a modest, though non negligible spatial 
polarization of the rotors is present also in the parame- 
ters region that was destined to homogeneous phases in 
the limiting case /i = 0. 

Starting from this setting, we have decided to perform 
a campaign of Vlasov based simulations to challenge the 
rich scenario predicted within the realm of the Lynden- 
Bell violent relaxation theory. By numerically solving the 
Vlasov equation, we avoid dealing with finite size eff'ects, 
as stemming in direct iV-body schemes, and so provide a 
more reliable assessment of the overall correctness of the 
theory. The results of the investigations are reported in 
the forthcoming sections. 
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III. MAGNETIZATION AND ITS 
FLUCTUATIONS 

The Vlasov equation ^ can be resolved numerically. 
To this end, we use the semi-Lagrangian method with 
cubic spline interpolation, as implemented in the vmf 90 
program that has been used already in Ref. with the 
HMF model. 

In order to study the properties of the QSS regime, we 
adopt the following procedure: 

1. The system is started with a waterbag initial con- 
dition (g]). 

2. It is run without collecting data between times to = 
and ti = 100. 

3. Time averages of the magnetization m, of rrix and of 
the variance of the magnetization m are performed 
in the time range between ti = 100 and t2 = 200, 
defining 



1 



t2-ti 



t2~h 



m{t)dt 



ix{t)dt , 



t2 



(m(t) - mfdt 



(8) 



(9) 



(10) 



We thus skip the strong oscillations of the initial vio- 
lent relaxation and focus on the subsequent dynamical 
regime, where, anyway, some oscillations are still present, 
which we quantify by the standard deviation am- We re- 
peat the above procedure on a grid of 39 by 39 points 
in the (too?"") plane, each one corresponding to a sim- 
ulation. Performing such a study allows us to assess in 
a systematic manner the behaviour of the average mag- 
netization in the QSS regime and to compare numerical 
results with Lyndcn-Bell's theory. Whether the theory is 
flexible enough to accommodate all the general features 
of the resulting diagram depends on it taking into ac- 
count in a comprehensive manner the behaviour of the 
model. 

The average value of the magnetization taken from 
Vlasov simulations is displayed in Fig. [2l The line of 
transition provided by Lynden-Bell's theory is displayed 
on top and we observe that it separates satisfactorily the 
region m > from the region m « for rriQ ^ O.6. The 
transition is sharp for low values of mg, corresponding 
to the prediction of Lynden-Bell's theory that the tran- 
sition is of first order. The transition is smoother for 
larger values of mg, corresponding to the Lynden-Bell's 
prediction of a second order transition. We display in 
Fig. [3] the same diagram as the one of Fig. [51 but us- 
ing rux instead of m. The general aspect of the diagram 




FIG. 2: Average magnetization fh for the HMF model with 
no external field. The transition predicted by Lynden-Bell's 
theory is indicated by the full black line. 



is similar, although, as remarked in [13| . the transition 
looks overall sharper when using rux as an order param- 
eter. We note that Lynden-Bell's transition line sepa- 
rates very well the non-homogeneous from the homoge- 
neous phase for uiq ^0.4. For higher values of mo there 
are simulations for which rrix < (blue spots below the 
Lyndcn-BcU's transition line in Fig.|3l). This occurs when 
the phase of m is tt (instead of zero). The phase could 
in principle take any value in [— tt, tt], but since my = 
the only dynamically accessible values are and tt. The 
fact that TTix can flip from positive to negative values is 
also an indication of the presence of a second-order phase 
transition. Indeed, these flips are not present in the first 
order phase transition region, where an cntropic barrier 
at the phase transition separates positive and negative 
values of nix ■ The difference in the values of m and nix 
can only arise from time fluctuations of ruxit). Indeed 

m? ~ rn^ ~ nix^ + cr^ (11) 

To illustrate the difference between these two quantities, 
we compare them in Fig. |4]for mo = 0.1 and 0.4. In the 
low energy phase the two quantities are indistinguishable, 
prooving that the fluctuations are small. It is conflrmed 
that nix goes sharply to zero at the transitions energy and 
remains zero in the whole high energy phase, as found in 
[l3| . On the contrary, m, the quantity measured in [§], 
has a tail of positive values at high energy, especially 
visible for mg = 0.4, proving that fluctuations are here 
larger. 

The variance of the magnetization, (7™, is displayed 
in Fig. [5l It is confirmed that, below the transition line 
predicted by Lynden-Bell's theory, fluctuations of m are 
small. They are instead large in the high energy region 
above the Lynden-Bell's transition line for toq > 0.4. 

By pointing out the different results that arise from 
the choice of different order parameters, m or nix, we 



5 




FIG. 3: Average of the a;-component of the magnetization nix 
for the HMF model with no external field. The transition pre- 
dicted by Lynden-Bell's theory is indicated by the full black 
line. 
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FIG. 4: Comparison of m and rrix as a function of u for 
mo = 0.1 and mo = 0.4. 

hope that in future studies the problem of the out-of- 
equilibrium phase transition in the HMF model will be 
analyzed more carefully. 



IV. RESPONSE TO THE APPLICATION OF A 
SMALL MAGNETIC FIELD 

In this Section, we present the results of simulations for 
the HMF model with a small external magnetic field, we 
choose h = 0.1. The average value of the magnetization 
TO obtained from Vlasov simulations is displayed in Fig. [51 
The phase transition is removed by the application of the 
field, as it happens for equilibrium phase transitions. The 



FIG. 5: Amplitude of oscillations measured by cr„i for the 
HMF model with no external field. 
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FIG. 6: Average magnetization m for the HMF model with 
an external field h = 0.1. 



magnetization, for all values of toq, decreases smoothly 
to zero as the energy is increased. 

Magnetic susceptibility x = dm/dh characterizes the 
response of the system to the application of an external 
field. It has been shown in Ref. [10] that certain param- 
eter regions display a negative magnetic susceptibility. 
This is a signature of ensemble inequivalence, shown here 
in a out- of- equilibrium setting, as the system is trapped 
in the QSS regime and docs not reach equilibrium. In 
this Section, we provide a similar measure by taking the 
difference of the average magnetization between simula- 
tions with h = 0.1 and simulations with h ~ 0. The 
result is displayed in Fig. [71 While our computations 
provide a discrete difference instead of a derivative, ob- 
taining a lower value of the average magnetization for 
h — 0.1 than for /i = is the sign of a negative suscep- 
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FIG. 7: Difference of average magnetization between simula- 
tions with h = 0.1 and simulations with h — 0. The white 
lines indicate the zero level, so that the darker region close to 
(mo = 0,u — 1/2) is a region of negative magnetic suscepti- 
bility. 

tibility nonetheless. As expected from Ref. a region 
of Fig. [7] displays x < for low values of mo, in the 
vicinity of the first-order transition found in the theory. 
We are thus able to confirm the theoretical prediction on 
the basis of Vlasov simulations. Figure [7] also displays a 
large value of x around the transition line predicted by 
Lynden-Bell's theory. 

V. CONCLUSIONS 

We have performed a study on the adequateness of 
Lynden-Bell's theory as compared to numerical simula- 



tions of the Vlasov equation for the Hamiltonian Mean- 
Field model. Our results confirm previous studies based 
on iV-body simulations on the general quality of the 
phase diagram. We extended the knowledge of the phase 
diagram by several additional measurements: the ampli- 
tude of oscillations, (T,„, the magnetic susceptibility. By 
doing so, we point out that in regions where non negligi- 
ble fiuctuations of the magnetization m occur, the theory 
is not expected to work, whereas the agreement is quite 
good between numerical simulations and theory in re- 
gions where fluctuations are small. We also confirmed 
that there are reg ions of negative susceptibility, as pre- 
dicted in Ref. jlOf . 

Finally, we also discuss the more fundamental, but re- 
lated issue, of the appropriate thermodynamical quan- 
tity to follow in the simulations. This latter issue is not 
touched upon in the literature but reveals qualitatively 
different results for the transition from magnetized to 
non-magnetized regimes. As of now, Lynden-Bell's the- 
ory provides a clear determination of the order of the 
transition and a tricritical point is found in the phase 
diagram. However, the dynamical aspects of the transi- 
tions are not yet elucidated, as is known from numerical 
simulations, even if steps are taken in that direction 
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